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1.  INTRODUCTION 
1.1  Background 

The  primary  objective  of  optical  devices  in  non-lethal  applications  is  to  elicit  a  temporary 
response  through  the  targeted  individual’s  visual  system.  While  low-level  exposures  are  useful 
for  providing  a  warning  response,  systems  that  produce  high-level  exposures  that  may  degrade 
visual  function  temporarily  are  increasingly  sought  after  for  tactical  advantages.  While  systems 
utilizing  either  broadband  optical  radiation  or  laser  sources  are  feasible  for  this  use,  functional 
application  of  both  technologies  and  quantifying  the  potential  degradation  are  key  steps  in  the 
weapons  development  process.  The  degradation  effects  require  high-fidelity  assessments  of 
exposures  on  human  vision.  To  meet  this  need,  the  Air  Force  Research  Laboratory,  711th 
Human  Performance  Wing,  Human  Effectiveness  Directorate,  Bioeffects  Division,  Optical 
Radiation  Bioeffects  Branch  (711  HPW/RHDO)  previously  developed  a  computational  model  to 
estimate  the  transient  visual  effects  created  by  exposure  to  high-intensity  light  sources.1 

The  current  model  was  an  evolution  of  a  conceptual  model  derived  largely  from  research  on  the 
human  visual  response  to  high-intensity  monochromatic  laser  light  sources  conducted  in  711 
HPW/RHDO  over  the  previous  decades.2,3  Those  data,  along  with  empirical  models  accepted  in 
the  scientific  community  for  determining  intra-ocular  light  scatter  and  contrast  sensitivity, 
resulted  in  a  first-order  model  of  glare  and  photostress  recovery.  The  original  model  was  image- 
based  and  took  a  radiometric-referenced  image  of  the  optical  source  (an  extended  source)  and 
processed  it  through  a  two-dimensional  intra-ocular  scatter  model  to  calculate  the  light 
distribution  on  the  retina,  superimposed  this  distribution  on  the  visual  task,  and  estimated  the 
glare  effect  induced  by  the  reduction  in  target  contrast. 

For  flash  recovery,  a  flashblindness  model  was  used  to  calculate  the  resultant  transient  relative 
scotoma  at  any  given  time  after  the  exposure;  this  scotoma  was  superimposed  on  the  visual  task 
to  estimate  the  reduction  in  the  observed  target  contrast.  For  both  glare  and  flash  recovery,  the 
visual  task  was  of  uniform  luminance  contrast,  and  obscuration  levels  were  calculated  by 
comparing  the  image  contrast  with  a  contrast  threshold  estimate,  the  contrast  threshold  being 
based  upon  the  size  of  the  visual  task  (angular  subtense)  and  the  ambient  luminance.1 
Obscuration  was  represented  visually  via  a  binary  illustration  of  the  portions  of  the  task  in  the 
scene  that  were  visible,  and  calculated  as  a  percentage  of  the  target  obscured.  Using  this 
approach,  predictions  of  relative  effectiveness  for  candidate  non-lethal  flashbang  devices  were 
performed  as  part  of  the  improved  flashbang  grenade  program.4,5 

Recently,  the  vision  model  was  enhanced  with  the  inclusion  of  an  approach  that  replicates  the 
spatial  frequency  processing  mechanisms  that  are  found  in  the  human  visual  cortex.6  The  model, 
called  PREMO  (for  Photostress  REcovery  MOdel),  approximates  cortical  processing  by  taking  a 

7  & 

complex  two-dimensional  visual  scene  and  decomposing  it  into  six  spatial  frequency  channels.  ’ 
This  bandpass-filtered  approach  is  founded  on  prior  research  published  in  the  vision  science 
literature  for  target  detection  and  identification.  ’  In  PREMO,  visual  recovery  following  a  flash 
exposure  is  estimated  by  modeling  recovery  for  each  spatial  frequency  channel,  and 
reconstructing  the  visual  scene  at  different  times  after  the  flash  by  recombining  results  from  the 
recovering  channels.  The  advantage  of  this  approach  over  the  prior  model  is  that  it  does  not 
compress  the  visual  world  into  a  single  size  and  contrast;  rather  it  represents  the  visual  task  in  all 
its  complexity,  comprised  of  a  multitude  of  target  “sizes”,  contrasts,  and  luminance  levels. 
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Image-based  sources  and  visual  tasks  supported  by  PREMO  allow  the  user  flexibility  in 
analyzing  a  wide  variety  of  exposure  conditions  and  sources.  For  the  source,  digital  images  from 
searchlights,  LED  arrays,  thermobaric  devices,  and  other  sources  can  be  used.  Because  PREMO 
applies  digital  imagery  as  the  visual  task,  tasks  are  no  longer  limited  to  simple  uniform  geometric 
shapes:  virtually  any  grayscale  digital  image  can  be  used.  In  addition,  the  results  are  presented 
as  a  probability  of  identification  at  any  given  time  after  exposure.  The  model  also  has  a 
visualization  window  that  renders  how  the  visual  image  is  obscured  by  the  optical  exposure  and 
how  vision  recovers  after  the  exposure  has  been  terminated.  These  features  provide  the 
capability  of  the  model  to  cover  a  wide  range  of  operational  scenarios  and  allow  the  user  to 
envisage  the  consequences  of  optical  radiation  visual  effects. 

1.2  Purpose  and  Scope 

The  purpose  of  this  report  is  to  document  the  mathematical  basis  of  the  Photostress  Recovery 
Model  (PREMO)  and  to  describe  its  logic,  function,  and  scope.  This  report  describes  the  model 
algorithms  in  detail  and  their  scientific  underpinnings.  Brief  appendices  provide  additional 
details  on  the  use  of  the  model  including  the  software  architecture  and  preprocessing  of  the 
optical  radiation  sources. 
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2.  MODEL  DESCRIPTION 
2.1  Overview 

PREMO  is  a  computational  model  that  predicts  the  impact  of  broadband  and  monochromatic 
optical  radiation  on  visual  function.  The  model  accepts  images  of  the  optical  source  and  the 
visual  task  to  analyze  the  glare  and  photostress  recovery  effects  (Figure  1).  In  practical  terms  the 
model  analyzes  how  visible  light  from  intense  optical  sources  affects  the  visibility  of  objects  in  a 
viewed  scene  and  the  time  course  of  image  detection  recovery  after  the  obstructive  optical  source 
has  been  removed. 


Flash  Image 


Visual  Task  Image 
Or  3D  Model 


Adjust... 

-  Field  of  View 

-  Position  of  Observer 

-  Position  of  Visual  Task  (VT) 

-  Orientation  of  VT  (model) 

-  Size  of  VT  (image) 

-  Position  of  Flash  Source 

-  Size  of  Flash  Source 

-  Position  of  Viewer 

-  Look  Vector 


Generate  two 
images 


Visible  Portions  of 
Visual  Task  at 
Time  (t) 


Figure  1.  An  illustration  of  how  the  model  combines  the  images  of  the  optical  source  and  the 
visual  task  to  analyze  the  glare  and  photostress  recovery  effects 

A  description  of  the  logical  steps  taken  by  the  model  follows: 

•  A  scene  generator  is  used  to  transform  the  physical  locations  of  the  optical  source, 
observer,  line  of  sight,  and  visual  task  into  a  two-dimensional  scene  from  the  observer’s 
point  of  view. 

•  The  model  input  parameters  include  the  source  parameters,  the  source  image,  the  mission 
scenario,  complex  scene  parameters,  a  three  dimensional  object  or  an  image  of  the 
complex  scene,  observer  parameters,  and  the  field  of  view.  For  an  extended  source  the 
input  parameters  include  the  size,  duration,  measured  target/object  illuminance,  and  range 
at  which  illuminance  was  measured.  For  a  laser  source  the  input  parameters  include  the 
wavelength,  beam  waist,  range  to  observer,  aperture  waist,  divergence,  and  laser  power. 
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For  the  complex  scene,  an  image  of  the  scene  is  used  as  input  and  the  user  specifies  the 
image  size  and  the  luminance  range  for  the  scene. 


•  The  glare  and  photostress  recovery  algorithms  process  the  source  input  parameters  to 
develop  a  retinal  luminance  image  distribution.  The  retinal  luminance  image  distribution 
can  be  produced  using  a  geometric  optics  approach  (recommended  for  extended  sources) 
or  convolution  of  the  corneal  luminance  with  the  point  spread  function  (PSF)  for  the 
human  eye.  The  resulting  distribution  is  processed  through  the  spatial  frequency 
algorithm. 

•  For  the  scene  image  the  model  generates  band-pass  luminance  distribution  matrices  and 
low-pass  luminance  distribution  matrices.  For  each  band-pass  luminance  distribution 
matrix  of  the  scene,  the  spatial  frequency  algorithm  computes  a  reduced-contrast  matrix 
and  compares  it  to  a  detection  threshold  to  determine  recovery. 

•  The  photostress  recovery  model  determines  recovery  of  the  scene  at  the  time  steps 
provided  by  the  user. 

•  Given  the  recovered  scene  at  a  given  point  in  time,  the  probability  of 
identification/recognition  algorithm  provides  the  probability  of  identifying/recognizing  a 
target  based  on  the  recovered  spatial  frequency  content  and  the  extent  of  detectable 
edges. 

A  detailed  description  of  the  mathematical  formulations  and  the  scientific  underpinnings  of  the 
model  are  presented  in  the  following  sections. 

2.2  Components 
2.2.1  Scene  generation. 

The  OpenSceneGraph'  3D  viewer  allows  the  user  to  position  one  or  more  optical  sources,  the 
observer,  and  the  visual  task  in  3D  space.  In  order  to  create  an  accurate  representation  of  the 
observer’s  field  of  view,  a  3D  virtual  representation  of  the  scenario  is  used  to  generate  input 
images  that  can  be  passed  and  processed  by  the  model.  This  is  illustrated  in  Figure  3.  This  3D 
virtual  representation  of  the  scene  requires  these  parameters: 

Source  Image.  This  image  file  is  a  square,  gray-scale  image  that  represents  the  spatial 
distribution  of  the  optical  source.  It  is  typically  created  as  a  two-dimensional  circular 
Gaussian  gradient  with  white  at  the  center  and  black  on  the  edges.  In  the  figure  the 
optical  source  is  represented  by  a  thermobaric  device  (flashbang).* 

Visual  Task  Image.  This  image  file  represents  the  visual  task  in  the  scene  that  is  to  be 
obscured.  A  3D  model  file  (e.g.  a  3D  model  of  a  soldier  created  in  3D  Studio  Max,  as  in 


1  See  Appendix  A  for  a  description  of  the  OpenSceneGraph  software  architecture  and  the  type  of  files  it 
supports. 

*  In  the  absence  of  a  user-supplied  source  image,  a  simulated  source  image  can  be  generated  by  the  model  (a 
Gaussian  or  a  flat  top  distribution  are  currently  available). 
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the  figure)  or  a  2D  image  may  be  used  to  represent  the  visual  task  in  the  3D  scene.  The 
vision  model  requires  luminance  values  for  each  pixel  in  the  complex  scene.  In  its 
current  configuration,  the  PREMO  accepts  maximum  and  minimum  luminance  values  as 
user  inputs,  and  uses  these  together  with  the  image  greyscale  values  to  scale  the 
luminance  of  each  pixel  within  the  complex  scene. 

Scenario  Geometry.  The  scenario  geometry  dictates  the  positioning  of  all  of  the  elements  in 
the  scene  and  includes 

a)  Visual  task  position,  orientation,  and  scale, 

b)  Source  image  position,  orientation,  and  scale, 

c)  Observer  position,  orientation,  and  field  of  view.  The  observer  is  shown  as  a 
schematized  camera  (in  yellow)  in  the  figure.  The  scene  from  the  observer’s  point  of 
view  is  shown  in  Figure  3. 


Figure  2.  A  visual  representation  of  the  scenario  in  3D  space 
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The  geometric  relationships  between  the  objects  in  the  scenario  determine  the  visual  effect  of  the 
source.  The  observer,  source(s),  and  visual  task  can  each  be  independently  placed  anywhere  in 
the  scene.  In  most  situations  the  source  needs  to  be  in  the  vicinity  of  the  line  of  sight  to  have  any 
significant  effect.  The  underlining  infrastructure  supports  a  different  point  of  fixation  for  each 
time  step  to  simulate  eye  movements,  but  typically  only  one  or  two  points  of  fixation  are  used  for 
the  entire  simulation. 


2.2.2  Optical  source  processing. 

This  component  calculates  the  optical  source  fluence  at  the  cornea  for  the  viewing  distance  and 
propagates  the  source  through  the  optical  media  of  the  eye  to  create  an  image  on  the  retina. 

2.2.2.1  Corneal  light  distribution. 

For  a  broadband  source,  light  from  the  source5  is  propagated  to  the  eye  to  obtain  the  corneal  light 
distribution.  The  quantity  source  can  either  be  a  single  value  (for  a  point  source  such  as  a  laser) 
or  an  array,  as  in  an  image  (for  an  extended  source).  To  calculate  the  corneal  light  distribution 
for  an  extended  source,  the  model  uses  an  image  of  the  source  and  its  illuminance,  E(D),  which 
is  an  input  value  to  the  model.  If  the  illuminance  was  originally  measured  at  a  given  distance,  D, 
the  inverse  square  law  is  used  to  evaluate  the  illuminance,  E(R),  at  the  observer  distance,  R: 


E(R)  =  E(D)  ■ 


D 

\Rj 


(1) 


This  computed  illuminance  is  then  used  with  the  source  image  to  define  the  corneal  illuminance 
distribution.  For  each  pixel  in  the  source  image  the  corneal  illuminance,  E(xi,yj),  (for  glare)  and 
the  integrated  illuminance,  H(xi,yjj ),  (for  photostress  recovery)  are  given  by  the  following 
equations: 


E(xl,yJ)  = 


EjRygjj^yj) 


H(xl,yJ)  = 


t-E(R)-g{xi,yj) 

2>(W;) 


(2) 

(3) 


where  g(xu  v  j  is  the  gray-scale  value  of  a  pixel  in  the  source  image. 

For  a  laser  source,  the  model  calculates  the  corneal  light  distribution  from  the  laser  parameters 
supplied  by  the  user;  wavelength,  divergence,  beam  waist,  range  to  the  observer,  beam  profile, 
and  distribution  (Gaussian  or  flat  top).  Given  the  laser  parameters,  the  beam  diameter  at  the 
observer  location  is  given  by  the  following  equation: 


§  The  specifications  of  the  optical  source,  including  the  required  input  parameters,  source  image,  and 
conversion  to  photometric  units  are  detailed  in  Appendix  B  for  both  laser  and  broadband  sources. 
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Dl  =  +  (r-r0)V2  • 


(4) 


The  parameter  a  is  the  beam  waist,  r  is  the  range  to  observer,  rQ  is  the  aperture  waist,  and  (j>  is 
the  beam  divergence. 

The  laser  power  transmitted  by  an  aperture  diameter,  D/,  is  given  by: 


(5) 

The  parameter  ®0  is  the  laser  output  power.  For  a  circular  beam,  the  irradiance  at  the  cornea  is 
given  by: 


d>  =  0 

o 


E  = 


0 


y 

j 


1.27Q 


The  Illuminance  Ev  is  related  to  the  spectral  irradiance  Ex  by  the  following  equation 


(6) 


Ev=  683  •  [  EA  -  VA  •  dA 

J  ,  (7) 

where  Vx  is  the  photopic  luminous  efficiency  function. 

2. 2. 2. 2  Intraocular  light  scatter. 

Light  from  a  point  source  is  not  focused  perfectly  on  the  retina,  but  instead  it  is  scattered  within 
the  eye,  which  produces  a  veiling  luminance  over  the  field  of  view.  The  scattered  light  causes  a 
reduction  in  image  contrast,  particularly  for  objects  in  the  field  of  view  close  to  the  source.  The 
masking  effect  caused  by  the  light  scattered  in  the  ocular  media  is  termed  veiling  glare.  To 
evaluate  the  effects  of  veiling  glare,  an  algorithm  for  intraocular  scatter  is  used  to  model  the 
distribution  of  a  light  source  on  the  retina.  The  function  used  to  predict  the  distribution  of  the 
light  scattered  based  on  the  light  distribution  at  the  cornea  is  referred  to  as  the  point  spread 
function  (PSF).10'13 

1  O 

An  analytical  expression  for  the  PSF  that  closely  fits  experimental  data  was  developed  by  Vos. 
The  Vos  formulation  distinguishes  the  different  components  of  entoptic  light  scatter  (anterior 
segment,  scleral  wall,  and  epithelium)  and  includes  the  dependency  of  the  PSF  on  age  (A)  and 
eye  pigmentation  (P).  The  current  model  uses  this  expression  to  compute  the  image  profile  of  a 
point  source  on  the  retina: 
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PSF{6 )  =  Anterior(9 )  +  Wall(9 )  +  Epithelium(9)  , 


(8) 


where: 
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(11) 


and  9  is  eccentricity  from  the  center  of  the  source  on  the  retina  in  degrees. 


Figure  4  shows  the  shape  of  the  PSF,  which  describes  the  magnitude  of  the  intraocular  scatter  as 
a  function  of  the  angle  between  the  incoming  source  and  the  retinal  location.  The  intensity  is 
greatest  for  the  central,  zero-degrees  location,  and  falls  off  as  the  angle  increases.  For  an 
extended  source  the  same  scattering  profile  is  used,  but  it  is  convolved  with  a  two-dimensional 
light  distribution  rather  than  a  point  source. 


Figure  4.  The  intensity  profile  of  a  point  light  source  as  a  function  of  retinal  angle  due  to 

optical  scatter  within  the  eye 
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The  model  convolves  the  input  comeal  light  distribution,  E(x,  y ),  with  the  intraocular  point 
spread  function,  PSF(O).  The  output  of  the  algorithm  is  a  retinal  image  luminance  profile: 

L{x,  y)  =  E(x,  y )  0  PSF{6 ) ,  (12) 

where  L(x,y)  is  the  retinal  luminance  distribution  in  cd-m’2,  which  represents  the  source  glare 
field. 

2.2.23  Geometrical  optics  approximation. 

For  an  extended  source,  light  from  the  flash  is  collected  by  the  eye  and  the  flash  image  is  focused 
onto  the  retina.  A  schematic  diagram  is  shown  in  Figure  5. 


Figure  5.  Light  collected  by  the  lens  and  focused  by  the  eye 

Light  with  luminance,  fry  is  emitted  from  an  element  in  the  flash  with  area  dAg.  The  distance 
from  the  element  to  the  eye  is  Rg.  Across  all  the  elements  in  the  flash,  the  average  luminance  of 
the  flash,  L,  for  N  elements  is  given  by: 


L 


(13) 


For  an  image  of  the  flash,  in  terms  of  gray  scale  gg,  the  average  gray  scale  value,  g,  is  given  by: 


g  = 


ly  hj 


(14) 


and  the  luminance  of  each  element  is: 


hj=L 


g 


«• J 


g 


(15) 


For  each  elemental  area,  dAg,  the  illuminance,  Eg,  at  the  eye  is  calculated  from  the  source 
measurements  scaled  for  the  viewing  distance  by  the  inverse  square  law: 
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The  total  illuminance  at  the  eye,  E,  is  known  from  sum  of  the  individual  elements, 


£  =  2X 

■  (17) 

Assuming  that  all  elements  have  the  same  elemental  area,  dA,  and  are  located  a  distance  R  from 
the  lens,  the  total  illuminance  is  given  by 


E  = 


NdA  A  T 
R2  ~  R2 


(18) 


Knowing  the  area  of  the  flash,  A,  and  the  distance  from  the  observer,  R,  the  retinal  luminance  of 
each  element,  Luj,  is  given  by 


L 


i-j 


?jE(R)S- 
A  g 


hj 


A  g 


(19) 


2.2.3  Spatial  frequency  hierarchy. 

Psychological  studies  have  demonstrated  the  existence  of  a  frequency  selective  mechanism  in  the 
human  visual  system.  Previously,  a  number  of  models  have  been  developed  to  simulate  the 
human  vision  system  response  using  a  set  of  cortex  filters.7,14.  Our  current  cortex  filters 
modeling  is  based  on  a  model  originally  developed  by  Watson14,  as  improved  by  Daly.15  As 
defined  by  both  Watson  and  Daly,  the  cortex  filters  are  modeled  by  two  classes  of  filters,  for 
radial  frequency  selectivity  and  for  orientation  selectivity.  Our  current  implementation  accounts 
for  the  radial  frequency  selectivity,  but  not  the  orientation  selectivity.  The  cortex  filters  are 
formed  as  differences  of  a  series  of  two-dimensional  mesa  filters  characterized  by  a  flat  pass- 
band,  a  transition  region,  and  a  flat  stop-band  region.  This  ensures  that  the  filter  set  sums  to 
unity.  The  transition  region  is  modeled  with  a  Hanning  window,  so  that  the  mesa  filter  can  be 
completely  described  by  its  half-amplitude  frequency ,pm,  and  a  transition  width,  tw,  as  follows: 


mesa  (p)  = 


1.0 


Vcos(^-^  +  <w/2)» 

2  tw 


0.0 


for  p  <  a 


tw 


/  2 


r  tw  tw 

for  A/2-y<P<A/2+y 


for  p  >  pU2  + 


tw 


(20) 


where  each  cortex  filter  is  characterized  by  a  flat  pass-band  and  a  transition  region.  The  cortex 
filter  is  formed  by  the  difference  of  two  mesa  filters  evaluated  with  different  half-amplitude 
frequencies.  The  kth  cortex  filter  is  given  by: 
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cortex k  (p)  =  < 


mesa(p)  | 
mesa(p)\ 


Pin  (^_1) 


Pm  (^-1) 


-mesa(p)\nAt) 
-base(p)  |A;j(t) 


for  k  =  l,  K -2 
for  k  —  K  —  \ 


(21) 


where  the  total  number  of  filters  is  K  and  the  half-amplitude  frequency  for  a  mesa  filter,  pi/2(k), 
is  given  by: 


P\n 


k 

(*)  =  2  . 


(22) 


The  lowest-frequency  fdter,  base(p),  referred  to  in  this  document  as  a  low-pass  filter,  is  modeled 
by  a  truncated  Gaussian  function,  given  by: 


e-{p2t  ICT1) 


base(p)  = 


0 


where  <r  is  defined  as: 


and  p//2  =  2~k. 


r  tw 

for  p<pl/2+— 

r  tW 

for  p>pl/2+—  t 


(23) 


(24) 


In  the  current  implementation,  the  transition  width,  tw,  of  each  filter  is  a  function  of  its  defined 
half-amplitude  frequency  as  follows: 


tw  =  -  ^1/2 . 

3 


(25) 


The  resulting  filter  set  is  shown  in  Figure  6. 


Figure  6.  The  cortex  filter  series 
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2.2.4  Glare  and  photostress  recovery. 

PREMO  predicts  two  effects  of  high- intensity  light  sources:  veiling  glare  and  photostress 
recovery.  Glare  is  defined  as  the  masking  effect  of  a  bright  light  caused  by  light  scattered  in  the 
ocular  media  which  produces  a  veiling  luminance  over  the  field  of  view.  Photostress  is  a 
psychophysical  process  which  persists  after  a  bright  light  exposure.  This  effect  is  because 
photostress  induces  an  afterimage  in  the  form  of  a  transient  scotoma  that  subsequently  recedes; 
the  effect  is  temporary  and  visual  sensitivity  eventually  returns  to  normal  and  normal  vision  is 
restored.  Photostress  recovery  time  (PSRT)  is  the  time  between  when  the  light  is  removed  and 
when  the  observer  can  see  a  particular  visual  stimulus.  In  the  study  of  both  glare  and  photostress 
recovery,  the  effect  is  measured  by  comparison  to  a  physical  light  source  that  produces  the  same 
visual  effect;  this  has  given  rise  to  the  generally  accepted  metric  as  the  equivalent  veiling  or 
equivalent  background  luminance  (EBL).11  In  addition,  both  glare  and  photostress  recovery 
effects  depend  not  only  on  the  parameters  of  the  inducing  light  source,  but  also  on  the  properties 
of  the  visual  task  and  the  visual  scene.  The  size  and  brightness  of  the  visual  task,  the  ambient 
illumination  level,  and  visual  task  contrast  threshold  all  play  a  critical  role  in  determining  the 
extent  of  glare  and  the  time  course  of  photostress  recovery. 

2.2.4.1  Contrast  threshold  estimation. 

Contrast  is  the  difference  in  luminance  and/or  color  that  makes  an  object  distinguishable.  In 
visual  perception,  contrast  is  determined  by  the  difference  in  the  color  and  luminance  of  the 
object  and  other  objects  within  the  same  field  of  view.  The  model  uses  luminance  contrast 
detection  as  the  critical  visual  function.  Luminance  contrast  detection  is  the  ability  to  perceive 
the  presence  of  two  objects  to  a  difference  in  luminance  between  the  two.  Luminance  contrast 
threshold  and  its  measurement  have  been  extensively  studied  in  the  vision  science  literature  and 
are  well  understood.16,17  Luminance  contrast  detection  is  also  consistent  with  the  analysis  of 
optical  radiation  visual  effects  because  the  effect  of  an  intense  light  source  is  to  reduce  the 
luminance  contrast  of  objects  and  therefore  reduce  a  subject  observer’s  visibility. 

1  O 

The  contrast  sensitivity  function  (CSF)  of  the  human  eye  developed  by  Barten  is  used  in 
PREMO  to  calculate  contrast  thresholds  based  on  the  scene  parameters  (contrast  threshold  is 
simply  the  inverse  of  the  CSF).  The  CSF  gives  the  contrast  sensitivity  as  a  function  of  spatial 
frequency.  Figure  7  shows  a  family  of  contrast  threshold  curves  calculated  for  a  field  of  view 
size  of  10°  x  10°  and  various  visual  task  luminance  levels.  As  shown  in  Figure  7,  the  contrast 
threshold  reaches  its  minimum  between  1  and  10  cycles  per  degree.  For  a  standard  observer  the 
CSF  is  given  by  Equation  (26)  below. 


CSF(f)  = 


5200  e 


-0.0016x /ftl+100/L) 


144 


63 


04  r  +  0.64x/  X-513-  + 


1 


2fo 


1  -  e 


-0.02  xfz 


(26) 


In  this  equation /is  the  spatial  frequency  in  cycles/degree,  L  is  the  luminance  in  cd-m’2,  and  Zo 
is  the  angular  object  area  in  square  degrees.  The  constant  5200  in  the  numerator  is  valid  for 
binocular  viewing.  For  monocular  viewing,  it  is  a  factor  V2  smaller,  or  3700. 


12 

Distribution  A:  Approved  for  public  release;  distribution  unlimited.  PA  Case  No:  TSRL-PA-15-0029 


10.000 


0.001  -I - , - , — 

0.1  1.0  10.0 

spatial  frequency  (cycles  per  degree) 

Figure  7.  Contrast  threshold  calculated  using  Barten’s  formula  for  a  field  size  of  10°  x  10° 


2. 2. 4. 2  Processing  glare  effects. 

As  described  earlier,  the  model  uses  the  visual  task  image  grayscale  levels  and  user  inputs  for  the 
maximum  and  minimum  luminance  to  calculate  luminance  values  across  the  visual  scene.  This 
step  renders  the  image  of  the  scene  into  units  appropriate  for  it  to  be  added  to  the  glare  retinal 
image  distribution.  The  superposition  of  glare  on  the  visual  task  is  accomplished  by  adding  the 
glare  retinal  image  luminance  distribution,  defined  by  Equation  (12),  or,  for  an  extended  source, 
by  Equation  (19),  to  the  visual  task  image  luminance.  The  combined  image  is  known  as  the  EBL 
reduced  visual  task  image  ERVT(x,y )  and  it  is  given  by: 

ER  VT  (x,  y )  =  EBL(x,  y)  +  VT  (x,  y) ,  (27) 

where  EBL(xy)  is  the  retinal  image  distribution  of  the  glare,  and  VT(xyt)  is  the  visual  task  image 
luminance.  ERVT(xy)  is  convolved  with  each  of  the  six  cortex  filters  described  above  to  produce 
a  set  of  five  band-pass  filtered  images  denoted  by  ak{xyh  and  a  low-pass  image  ao(xy).  Each  band¬ 
pass  filtered  image,  ak(x<y) ,  is  obtained  using  Equation  (28): 

ak(x,y)  =  ERVT(x,y)*gk(x,y),  (28) 

where  the  array  gk(x,y)  represents  the  inverse  Fourier  Transform  of  the  cortex  filter,  cortex,,  (p) , 
and*  is  the  convolution  operator. 

The  low-pass  image,  a0(x,y),  is  generated  by  convolving  the  ERVT  image  with  the  lowest- 
frequency  filter  base(p), 


a0  (x,  v)  =  ER  VT (x,  y)  *  base(x,  v) 

? 
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(29) 


where  base(x,y)  is  the  inverse  Fourier  Transform  of  the  cortex  filter  base(p).  For  each  band-pass 
filtered  image  we  define  a  corresponding  local  luminance  mean  image,  l(x,y),  which  is  a  low- 
pass  filtered  version  of  the  image  containing  all  the  energy  below  the  band: 


i- 1 

Ii(x,y)  =  a0(x,  y)  +  D 


k= 1 


(30) 


A  two  dimensional  local  contrast  array,  Ci(x,y),  is  also  defined  for  each  band-pass  filtered  image. 
Contrast  is  computed  as  the  ratio  of  band-pass  filtered  image  to  the  corresponding  low-pass 
filtered  image: 


Ci(x,y) 


<*,-(*>  D 

l,(x,y) 

9 


(31) 


where  k(x,y)  >  0. 

This  definition  provides  a  local  contrast  measure  for  every  band  that  depends  on  the  local  energy 
at  that  band  and  the  local  background  luminance.  The  contrast  at  each  pixel  is  compared  to  the 
corresponding  contrast  threshold.  The  model  transforms  below-threshold  pixels  to  average 
luminance  to  show  which  part  of  the  band-pass  filtered  image  is  being  obscured  by  the  glare 
source.  The  recovered  scene  is  computed  by  adding  the  five  band-pass  filtered  images  each  five- 
threshold  transformed  and  the  low-pass  image  a0(x,y) . 

2. 2. 4. 3  Processing  ph  otostress  recovery  effects. 

The  photostress  recovery  analysis  describes  the  recovery  of  visual  sensitivity  after  the  source  has 
been  turned  off  though  a  visual  process  known  as  dark  adaptation.  As  with  the  glare  model,  the 
photostress  recovery  model  also  uses  the  equivalent  background  principal  to  model  the  recovery 
in  sensitivity.  The  effects  of  the  source  are  again  represented  as  a  veil  of  light,  but  the  effect 
fades  with  time.  The  initial  step  in  calculation  involves  expressing  the  flash  exposure  in  terms  of 
the  integrated  retinal  illuminance.  For  this,  by  analogy  to  the  glare  model,  the  PSF  is  applied  to 
the  flash  energy  to  estimate  the  retinal  distribution  of  the  flash,  but  calculated  as  an  integrated 
luminance  in  cd  m'2-s;  the  product  of  the  equivalent  background  luminance  and  the  exposure 
duration.  Then  the  integrated  luminance  is  expressed  in  terms  of  integrated  retinal  illuminance, 
RI,  in  photopic  troland-seconds  (td-s)  by  multiplying  by  the  area  of  the  pupil: 

RI(x,y)  =  L(x,y)  •  t  •  PA ,  (32) 

where  L(x,  y)  is  the  retinal  luminance  distribution  in  cd-m’2,  PA  is  the  pupil  area  in  mm2,  and  t  is 
the  exposure  duration  in  seconds.  RI(x,y)  is  the  resultant  integrated  retinal  illuminance 
distribution  for  the  exposure.  The  photostress  recovery  analysis  uses  the  integrated  retinal 
illuminance  from  the  photostress  stimulus  to  model  the  recovery  of  visual  sensitivity  over  time.19 
Recovery  is  represented  as  an  equivalent  background  luminance  distribution  that  decays  with 
time  according  to  a  power  function,  EBL(x,y,t).20  For  the  recovery  function,  the  integrated 
retinal  illuminance  distribution  is  first  converted  to  a  “flash  energy”  distribution,  F(x,y),  by 
taking  its  logarithm: 
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F(x,y)  =  log10[i?/(x,j)]. 


(33) 


These  values  are  then  used  as  input  into  the  decay  function,  and  at  any  time  after  exposure,  the 
EBL(x,  y,  t)  of  a  pixel  in  the  retinal  illuminance  distribution  is  given  by: 

EBL(x,  y,t)  =  1  Om'Loglo(t)  +  c ,  (34) 

where  m  and  c  are  functions  of  the  flash  energy  and  are  given  by: 

fl.OlF-  8.919  forF<  6.11  logtd-s 
m  =  <  (35) 

[-2.73224  forF>  6.11  logtd-s 

c  =  -0.4066-F2  +  6.4012  F- 19.427.  (36) 

Note  that  the  decay  function  takes  two  different  forms,  one  for  exposure  energy  values  less  than 
6.11  log  td-  s,  and  another  for  values  greater  than  this. 

The  behavior  of  the  function  is  shown  in  Figure  8,  which  shows  the  decay  of  the  EBL  with  time 
for  different  energies  in  the  initial  flash  exposure. 


Figure  8.  The  decay  of  the  EBL  after  the  exposure  is  terminated 

Recovery  of  visual  function  is  processed  differently  compared  to  the  glare  analysis.  This  is 
because  the  persistent  effect  of  the  flash  varies  with  time  after  the  initial  exposure.  For  the  glare 
the  EBL  is  used  as  a  static  contrast-reducing  mechanism,  whereas  for  photostress  recovery  the 
effect  of  the  EBL  at  time  t  must  be  evaluated  as  time  after  the  exposure  progresses. 
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For  the  recovery  analysis,  the  visual  task  image,  VT(x,y),  in  cd-m'2  is  convolved  with  each  of  the 
six  cortex  filters  described  previously  to  produce  a  set  of  five  band-pass  filtered  images  denoted 
by  cik(x,y)  and  a  low-pass  image  ao(x,y).  The  band-pass  filtered  image  ak(x,y,),  resulting  from  the 
convolution  is  written  as: 


ak  (x,y)  =  VT(x,y)*gk(x,y)  ^  ^ 

where  the  array  gk(x,y)  represents  the  inverse  Fourier  Transform  of  the  cortex  filter,  cortexfip), 
and  *  is  the  convolution  operator.  The  low-pass  image  ao(x,y )  is  generated  by  convolving  the 
visual  task  image  with  the  lowest-frequency  filter  base(p)\ 


«o  (x>  y )  =  VT (x’  y)  *  base(x,  y ) 


(38) 


where  base(x,y)  is  the  inverse  Fourier  Transform  of  the  cortex  filter  base(p). 

For  each  band-pass  filtered  image  a  two-dimensional  contrast  array  Ci(x,y,t)  is  defined.  For  each 
band-pass  filtered  visual  task  image  and  for  the  low-pass  visual  task  image,  the  contrast  is  given 
by  Equation  (38). 


C,(x,y,t) 


_ a,0  ,y) _ 

LPVTt(x,y )  +  EBL(x,y,t) 


i  =  0, 1, 2, 3, 4,  and  5 ,  (39) 


where,  LVPTi(x,y)  and  EBL(x,y,t)  denote  the  low-pass  visual  task  luminance  and  the  EBL  at  time 
t,  respectively. 

The  low-pass  visual  task  luminance  (LVPTfix.y))  for  a  given  visual  task  band-pass  image  is 
defined  as  the  visual  task  image  containing  all  energy  below  that  spatial  frequency  band. 


LPVTt(x,y ) 


i- 1 


a0(x,y)  +  Yjak(x,y)  i  =  2, 3, 4, 5 


k= 1 


a0(x,y) 


i  =  1 


(40) 


The  contrast  at  each  pixel  is  compared  to  the  corresponding  contrast  threshold.  The  model 
transforms  below-threshold  pixels  to  average  luminance  to  show  which  part  of  the  band-pass 
filtered  image  is  not  visible.  The  recovered  scene  is  computed  by  adding  the  six  band-pass 
filtered  images  with  individual  thresholds.  The  model  can  process  photostress  recovery  over  a 
period  of  time  to  show  a  time-elapsed  sequence  of  recovery.  Figure  9  illustrates  this  process. 
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Time  -  21  seconds  Time  -  31  seconds  Original  Image 

Figure  9.  An  illustration  of  how  recovery  of  visual  sensitivity  after  exposure  to  a  bright  light 

flash  is  represented  in  the  model 


Using  the  example  illustrated  in  Figure  9,  careful  inspection  of  the  recovery  images  reveals  a 
systematic  recovery  process  that  starts  with  low  spatial  frequencies  and  progresses  through  to  the 
high  frequencies  with  time.  The  scene  appeared  to  have  fully  recovered  to  include  most  of  the 
high  spatial  frequency  content  (fine  details)  after  3 1  seconds.  This  exemplifies  how  the  recovery 
time  for  a  given  visual  task  in  a  given  scene  will  be  highly  dependent  on  the  spatial  frequency 
content  of  the  scene.  Using  the  same  input  parameters,  Figure  10  shows  the  recovery  images  for 
each  of  the  spatial  frequency  bands. 
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Figure  10.  An  illustration  of  the  recovery  process  for  each  frequency  band 
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2.2.4.4  Visibility  metric. 

Developing  a  visibility  metric  for  arbitrary  task  located  in  an  arbitrary  scene  is  a  very 
challenging  task.  This  is  due  to  the  recovery  dependence  on  the  spatial  frequency  content  of  the 
scene.  An  initial  approach  was  taken  for  the  model  that  involved  assigning  a  visibility  weight  for 
each  filtered  band  and  assigning  a  contribution  share  as  a  fraction  for  each  band.  The 
obscuration  level  at  a  given  post-exposure  time  was  determined  from  the  number  of  pixels 
recovered,  the  weight  assigned  to  the  frequency  band,  and  the  value  of  the  contribution  share. 

2.2.5  Probability  of  identification/recognition. 

The  Visual  Effects  model  provides  two  approaches  to  predict  the  probability  of  identifying  a 
target  in  the  scene  at  any  given  post-exposure  time:  the  empirical  Target  Transform  Probability 
Function  (TTPF)  approach  developed  by  Johnson  and  Fawson21,  and  the  Contrast  Based  Edge 
Detection  approach  based  on  the  work  of  Johnson.22 

2.2. 5.1  Target  transform  probability  function  (TTPF). 

The  TTPF  relates  the  probability  of  performing  a  task  to  the  ratio  of  the  number  of  cycles  one 
can  resolve  across  the  minimum  dimension  of  a  target  to  the  standard  number  given  in  the 
Johnson  Criteria.23  Therefore,  given  the  level  of  resolvable  spatial  frequency  on  the  target  (N) 
and  the  recommended  spatial  frequency  level  for  50%  probability  to  detect,  identify,  or 
recognize  a  target  (N5o),  the  probability  can  be  predicted  based  upon  the  following  empirically 
derived  formulation: 


P,  = 


where 


(41) 


E 


2.7  + 0.7  x 


(  N  A 


\N50J 


(42) 


2.2.5. 1  Contrast  based  edge  detection  approach. 

The  formulation  of  an  edge  detection  algorithm  to  respond  to  contrast  rather  than  absolute 
luminance  values  provides  a  consistent  ground  for  edge  detection  in  a  scene  with  irregular  visible 
spectrum  distribution.  Given  the  luminance  distribution  of  the  scene,  the  contrast  based  Sobel 
operator  is  given  by  the  following  formula: 
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where  Ny  is  the  numerical  luminance  value  of  the  row  j,  column  i  element  of  the  matrix  that 
represent  the  scene.  The  symbol  ®  represents  the  convolution  operator,  while  the  convolution 
masks  denoted  by  A3  and  B3  are  given  by: 
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and  the  averaging  mask  C?  is  given  by: 

fl  1  A 

C3  =  1  1  1 


(44) 


(45) 


The  value  of  d  is  to  be  chosen  to  provide  an  acceptable  trade-off  between  noise  sensitivity  and 
the  dynamic  range  for  the  combination  of  the  human  visual  system  and  scene.  In  the  current 
implementation  the  value  of  d  is  set  to  unity,  but  a  threshold  formula  (Sth)  for  the  contrast-based 
Sobel  operator  as  a  function  of  contrast  (C)  was  developed  for  scene  luminance  ranges  from  2  to 
180  cd  m'2  (Equation  46  and  Figure  11).  Further  work  is  needed  to  extend  this  function  and 
define  thresholds  for  very  bright  and  very  dim  conditions. 


Slh  =  1 .5222  •  C4  - 3.9496  •  C3  +  2.959 1  •  C2  - 0. 1 565  •  C  +  0.3 1 7 


(46) 


Edge  Detection  Threshold  vs.  Contrast 


Figure  11.  Sobel  operator  threshold  as  a  function  of  contrast 

The  probability  of  finding  a  task  in  a  scene  at  a  given  post-exposure  time  can  be  predicted  using 
the  cumulative  normal  distribution,  cdf  (Equation  47),  given  that  the  fraction  of  recovered  edges, 
/i,  which  produces  recognition  in  50%  of  the  trials  is  established: 
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3.  APPLICATIONS 

This  experimental  approach  and  computational  analysis  model  can  be  used  for  a  variety  of 
applications  such  as  aviation  warning  systems,  military  operations,  and  law  enforcement  where 
intentional  optical  exposures  are  used  as  a  non-lethal  deterrent.  For  example,  visual 
effectiveness  of  non-lethal  dazzlers, 24-26  thermobaric  devices,5  and  laser  countermeasures  in 

97  98 

military  operations  have  been  assessed  with  the  more  simplistic  predecessor  of  this  model.  ’ 
The  model  may  also  be  used  to  evaluate  risks  to  safety  and  job  performance  in  situations  where 
intense  optical  sources  are  employed. 

4.  SUMMARY 

The  Photostress  Recovery  Model  (PREMO)  is  a  computational  model  of  the  glare  and 
photostress  recovery  effects  of  optical  radiation.  The  main  parts  of  the  model  are  a  scene 
generation  component,  an  optical  source  processing  component,  a  spatial  frequency  processing 
component,  a  glare  and  photostress  recovery  processing  component,  and  a  probability  of 
identification/recognition  processing  component.  The  optical-source  component  processes  input 
from  broadband  or  laser  sources  and  propagates  light  from  the  source  to  the  retina.  A  scene 
generation  module  specifies  the  spatial  relationships  between  the  observer,  radiation  source,  and 
visual  task  and  creates  a  representation  of  the  scene  from  the  observer’s  point  of  view.  The 
spatial  frequency  processing  component  combined  with  the  visual  effects  processing  component 
estimates  the  effects  of  glare  and  photostress  recovery  on  task  visibility  considering  the  spatial 
frequency  content  of  the  scene.  The  probability  of  identification/recognition  processing  provides 
an  estimate  of  the  probability  to  identify/recognize  a  target  in  the  scene  at  a  given  point  in  time. 
The  model  output  provides  a  quantitative  analysis  of  visual  obscuration  and  recovery.  The 
model  is  based  on  an  extensive  literature  on  optics,  vision  science,  and  digital  image  processing. 

PREMO  is  a  useful  tool  to  assess  the  extent  of  visual  disruption  for  a  variety  of  sources  and 
operational  scenarios.  The  current  instance  of  the  model  moves  beyond  simple  contrast  detection 
using  a  single  background  and  visual  task  luminance  into  contrast  detection  of  a  complex  scene 
with  varying  luminance  based  on  a  model  of  the  spatial  frequency  processing  in  the  visual 
cortex.  It  establishes  a  firm  foundation  for  future  research  for  providing  a  higher-fidelity  model 
and  understanding  of  the  temporal  effects  of  optical  radiation  on  the  human  vision  system. 

Separate  publications  that  describe  the  validation  studies  that  were  conducted  in  conjunction  with 
the  development  of  the  model  have  been  issued.29-31 
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Appendix  A:  Software  Architecture 


Flashblindness,  dll 

This  dynamically-linked,  Win32  library  written  in  ANSI  C++  currently  encapsulates  the 
computational  engine  for  computing  the  flashblindness  (transient)  effects  of  an  optical  source. 
This  library  depends  on  NumericalRecipes.dll. 

NumericalRecipes.  dll 

NumericalRecipes.dll  contains  code  based  on  the  book  "Numerical  Recipes  in  C++".  This  code 
contains  the  convolution  code  required  in  the  intraocular  light  scatter  model.  Flashblindness.dll 
also  uses  some  of  the  base  data  types  from  NumericalRecipes.dll. 

Open  Scene  Graph 

The  3D  functionality  of  the  Vision  Effects  library  is  provided  by  an  Open-Source  Scene  Graph 
library  called  OpenSceneGraph  (OSG).  OSG  is  an  objected  oriented  interface  to  the  OpenGL 
application  programming  interface  that  provides  the  3D  application  programmer  an  Object 
Oriented  interface  that  also  provides  most  of  the  necessary  infrastructure  required  for  a  3D  visual 
simulation.  OSG  makes  the  3D  environment  more  intuitive  and  object  oriented  and  makes  the 
scenes  more  efficient  by  limiting  costly  context  switching  in  OpenGL.  OSG  also  allows  3D 
model  and  image  files  of  various  formats  to  be  read.  The  types  of  files  that  can  be  read  are  based 
on  the  available  plug-ins.  These  plug-ins  have  the  format  osgDB_*.dll.  Looking  through  the  bin 
directory  will  reveal  the  types  of  models  and  images  that  can  be  loaded.  Most  importantly,  this 
plug-in  can  load  3D  models  of  the  format  *.3ds,  *.lwo,  *.lws,  *.dxf,  *.flt,  and  *.osg.  It  can  load 
image  files  of  the  type  *.bmp,  *.jpeg,  *.tiff,  and  *.png. 

Client  Applications 

A  plug-in  for  the  AFRL/RHDO  Application  Framework,  called  Flashblindness_Plugin.dll  was 
created  to  provide  a  more  consistent  and  feature  rich  user  interface.  As  the  framework  is  a  .NET 
application,  the  plug-in  assembly/DLL  was  written  in  C#.  An  intermediate  bridge  library, 
Flashblindness_Bridge.dll  was  written  in  managed  C++  to  bridge  managed  and  un-managed 
code. 
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Appendix  B:  Optical  Radiation  Sources 

Optical  sources  can  be  either  broadband  or  laser  devices.  Laser  beam  fluence  is  normally 
described  in  terms  of  W  m'  for  CW  exposures  or  J  m'  for  pulsed  exposures.  However,  for 
visual  effects  analysis,  laser  emission  must  be  converted  from  radiometric  to  photometric  units  to 
account  for  the  differential  spectral  sensitivity  of  the  eye.  The  spectral  luminosity  function  for 
photopic  vision  is  shown  in  Figure  11.  If  laser  fluence  is  specified  by  the  user  in  radiometric 
units,  the  model  converts  to  photometric  units  using  the  following  method. 


Wavelength  (nm) 

Figure  A-l.  The  CIE  photopic  spectral  luminosity  curve.  The  curve  shows  that  the  strength 
of  the  response  of  the  normal  human  eye  depends  on  the  incident  wavelength.  The  eye’s 
response  is  strongest  to  555  nm  under  photopic  or  daylight  conditions. 


The  photometric  equivalent  of  irradiance  is  illuminance,  Ev  (in  lm-m"2  =  lux),  which  is  related  to 
irradiance,  Ee  (in  W  m"2),  by: 


Ev=6S3xEexVA 

where  683  is  the  luminous  efficacy  constant  in  units  of  lm-W'1,  and  V\  is  the  photopic  spectral 
luminosity  coefficient.  In  addition,  since  the  human  visual  system  perceives  a  flickering  light 
source  with  a  frequency  above  60  Hz  as  continuous,  the  effects  of  pulsed  laser  systems  with 
pulse  repetition  frequencies  greater  than  60  Hz  can  also  be  analyzed.  For  these  pulsed  sources, 
the  irradiance  is  given  by  the  time-averaged  irradiance. 

The  above  relation  is  the  photometric  conversion  for  CW  sources  for  glare  analysis.  The  flash 
recovery  model  is  used  to  analyze  the  recovery  of  visual  sensitivity  after  the  source  has  been 
turned  off.  This  model  is  generally  appropriate  for  pulsed  or  short  duration  laser  exposures. 
However,  as  with  glare  sources,  multiple-pulse  exposures  greater  than  60  Hz  may  also  be 
analyzed  as  a  single  continuous  pulse.  For  pulsed  sources  the  relation  involves  the  conversion 
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from  radiant  energy  to  time-integrated  illuminance.  The  model  accepts  pulsed  laser  input  in 
terms  of  a  radiant  exposure,  He  (J-m~2).  The  radiant  exposure  is  converted  to  photometric  units 
by  the  equation: 


Hv=6S3xHexV^  (2) 

where  Hv  is  the  integrated  illuminance  in  lm-m  2-s. 

Broadband  optical  sources  are  processed  through  a  similar  procedure  except  that  the  spectral 
and,  if  necessary,  the  spatial  distribution  of  the  source  must  be  considered. 

The  eye  is  sensitive  to  optical  radiation  over  a  bandwidth  of  400-700  nm.  Therefore,  the  total 
source  irradiance  has  to  be  calculated  over  this  band.  The  model  accepts  a  text  file  of  the  relative 
spectral  distribution  of  the  source.  The  text  file  has  a  two  column  format  of  wavelength  number 
and  relative  intensity.  The  model  truncates  the  distribution  to  the  visible  band  and  interpolates 
between  wavelengths,  if  necessary.  The  user  may  also  choose  an  appropriate  spectral 
distribution  based  on  blackbody  color  temperature  from  a  built-in  database  if  no  spectral  data  are 
available. 

The  model  uses  the  source  or  irradiance  to  scale  the  relative  spectral  distribution  to  spectral 
irradiance  in  terms  of  W-m’2  ^'1.  The  spectral  irradiance,  Ex,  is  weighted  by  the  photopic 
luminosity  coefficient  at  each  wavelength,  and  is  integrated  across  the  visible  band  to  obtain 
illuminance,  Ey>  by: 


Ev  =  683  x  |  Ex  x  Vz  x  dA , 


(3) 


where  683  is  the  luminous  efficacy  constant  in  units  of  lm-W'1  and  V\  is  the  photopic  luminosity 
coefficient.  The  source  illuminance  Ev  is  then  multiplied  by  the  exposure  duration  t  to  obtain  the 
integrated  luminous  exposure  Hv  in  lm-m‘2-s_1  for  pulsed  sources 


H=E-t 


(4) 


If  the  source  irradiance  was  measured  at  a  distance  other  than  the  observation  distance  then  the 
luminous  exposure  has  to  be  corrected  for  viewing  distance.  For  a  source  measured  at  a  distance 
d,  the  inverse-square  law  is  used  to  calculate  the  luminous  exposure  at  the  observation  distance 
do/ and  is  given  by 


( 


V 


(5) 


Thus,  the  luminous  exposure  represents  the  total  brightness  of  the  source  at  the  observation 
distance.  In  the  case  where  the  source  is  extended,  that  is  occupies  space,  the  luminous  exposure 
must  also  be  scaled  by  the  physical  area  of  the  source.  The  physical  area  is  defined  by  the  user 
in  terms  of  a  length  and  a  width.  For  uniform  sources  the  luminous  exposure  is  simply  applied  to 
the  entire  surface.  However,  for  a  non-uniform  source  the  individual  elements  of  the  surface 
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must  be  calibrated  by  the  total  luminous  exposure.  This  is  accomplished  by  employing  the  user- 
supplied  2D  spatial  distribution  or  a  pre-defmed  Gaussian  spatial  distribution  supplied  by  the 
model.  This  distribution  has  parameters  for  both  the  horizontal  and  vertical  directions.  In  one 
dimension,  the  distribution  g(x)  has  the  form: 

g(x)  =  ex ;pf-(*-*o)  /A  (6) 

l  /(J  ) 

where  xo  and  crare  the  center  and  the  1/e  location  of  the  distribution,  respectively.  Whether  the 
flash  distribution  is  supplied  by  the  user  or  the  model,  the  calibration  procedure  is  the  same.  For 
an  m  x  n  array,  where  each  element  i  represents  the  relative  source  intensity  at  that  point  in 
space,  i  is  scaled  by  the  ratio  of  the  total  luminous  exposure  to  the  sum  of  the  relative  intensities 
in  the  array  to  obtain  the  calibrated  luminous  exposure  Hm  n 

II  _  i mj!  ^ v 

m,n  ~  m,n 

li 

1  •  (V) 
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